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Since it was theorized by Kerr in 1963, determining the spin of black holes from 
observed data was paid very little attention until few years back. The main reasons 
behind this were the unavailability of adequate data and the lack of appropriate 
techniques. In this article, we explore determining/predicting the spin of several 
black holes in X-ray binaries and in the center of galaxies, using X-ray and gamma- 
ray satellite data. For X-ray binaries, in order to explain observed quasi-periodic 
oscillations, our model predicts the spin parameter of underlying black holes. On 
the other hand, the nature of spin parameters of black holes in BL Lacs and Flat 
Spectrum Radio Quasars is predicted by studying the total luminosities of systems 
based on Fermi 7-ray data. All sources considered here exhibit characteristics of 
spinning black holes, which verifies natural existence of the Kerr metric. 
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I. INTRODUCTION 



The solution of Einstein equation for gravitational field of a spinning mass was obtained 



by Kerr which was later generalized by Boyer & Lindquist 
analytic extension of it. In the metric, the spin parameter 'a' of the gravitational object, 
i.e. black hole, is defined in such a way that |a| < M, when M is the mass of the black hole 
(speed of light (c) and Newton's gravitational constant (G) are chosen to be unity). Note 
that in the Boyer-Lindquist coordinate, the outer radius of the rotating, uncharged black 



2| to construct a maximal 



hole is Tj. = M + a/M2 - a2. Hence, for | a\ > M, the collapsed object will form a naked 
singularity, rather than a black hole, without an event horizon. In addition, for a = 0, the 
solution reduces to the Schwarzschild solution describing non-rotating black holes. Hence, 
predicting 'a' of black holes from observed data enables one to verify natural existence of 
the solution of Einstein equation for gravitational field of a spinning mass, more precisely 
the Kerr metric. 

Most of the black holes are either in X-ray binaries (XRBs) or at the center of galaxies. 
In XRBs they have masses at the most few tens of that of Sun, hence they are called stellar 
mass black hole (StBH). Examples of such black holes are GRS 1915+105 and GRO J165540. 
The black holes at the centre of galaxies have masses atleast million times that of Sun and 
are called super-massive black holes (SuBH). The black hole Sgr A*, at the centre of our 
galaxy, is an example of a SuBH. 

The aim of the present paper is to predict the spin of both kinds of above mentioned black 
holes. The value of a is expected to vary across the black holes, which is presumably being 
reflected in their various properties, e.g. quasi-periodic oscillations (QPOs), outfiows/jets, 
spectra, luminosities etc. Currently, the major methods of measuring a for black holes in 
XRBs are: (1) continuum spectrum fitting ^], (2) determining the shape of the gravitation- 
ally redshifted wing of an iron line 14], and (3) modeling QPO of the source which is supposed 
to be linked with the black hole [5|. For the present purpose, we adopt the third method 
in order to determine a in XRBs. As the QPO frequencies seem to be scaling inversely as 
mass (M) of the black hole, for SuBHs they are expected to be several orders of magnitude 
smaller than those in XRBs [G] and hence it is difficult to identify them in observed data. 
As a result, predicting a is uncertain for SuBHs based on QPOs. Therefore, we attempt to 
predict the spin of SuBHs, for the present purpose, based on underlying jet properties. In 
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particular, we adopt a class of AGNs with jets nearly aligned to the observer's line of sight, 
viz. blazar. 

In the next section, first we outline the properties of QPOs, then we model to describe 
them and finally predict spin of the black holes in XRBs. In §3, we briefly discuss the 
properties of blazars and then argue that the spin difference between Flat Spectrum Radio 
Quasars (FSRQs) and BL Lacs is responsible for the difference in their luminosities. Based 
on that we attempt to predict possible values of spin of black holes in FSRQs and BL Lacs. 
Note that blazars are radio loud active galaxies having very high luminosities which are 
expected to harbor a highly spinning black hole generally {3, 0] • Finally we summarize our 
results in 54. 



II. SPIN OF BLACK HOLES IN X-RAY BINARIES 



A. Classification of QPOs 

We attempt to predict the spin of the black holes in XRBs addressing the underlying 
QPOs. The QPOs in XRBs (particularly low mass XRBs) are generally of three types in 
terms of their observed frequencies — the ones of the order of a few hundreds to kilo Hertz: 
often called high frequency QPOs (HFQPOs) for black holes and kilo Hertz (kHz) QPOs for 
neutron stars, the others of the order of a few tens of Hertz: intermediate frequency QPOs 
(IFQPOs), and finally the QPOs of the order of one tenth of Hertz to few tens of Hertz: low 
frequency QPOs (LFQPOs). 



The HFQPOs/k 
pact object itself 



3z QPOs are expected to be linked with the spin parameter of the com- 
, |9[. In certain black hole and neutron star systems they appear in pair, 
e.g. GRS 1915+105, XTE J1550564, GRO J165540, 4U 1636-53, 4U 172834, Sco X-1. It was 
observed that frequencies in kHz pairs are separated by the order of either the spin frequency 
or half the spin frequency of the neutron star. Moreover, HFQPOs in pairs appear in a 2 : 3 



ratio. L 
as well 



owever, sometimes kHz QPOs, e.g. that in Sco X-1, seem to appear in a 2 : 3 ratio 



10| . All these properties favor the idea of a resonance mechanism behind the origin 
of HFQPOs/kHz QPOs. The difference between HFQPOs in black holes and kHz QPOs in 
neutron stars might be due to the presence of a rotating magnetosphere in the latter case, 
imprinting directly the spin frequency on the oscillations of the disk. In order to model their 



properties, earlier Mukhopadhyay js, 6| predicted spin of those underlying black holes and 
neutron stars, whenever unknown, for a given mass of the compact object. 

On the other hand, LFQPOs are thought to be produced by shock oscillations in accretion 
flows Uj and, unlike HFQPOs, expected not to be related directly with the spin of compact 
objects. Of course, with the change of spin, properties of accretion flows and consequently 
the underlying shock change, which can influence LFQPOs indirectly. What about IFQPOs? 
We will discuss them based on our model. 



B. Modeling QPOs 



Observations [12| indicate a strong correlation among the QPO frequencies in a wide 
range, which supports the idea that QPOs are universal physical processes. It was already 
shown that QPOs may arise from the nonlinear hydrodynamic/hydromagnetic phenomena 
of accretion flows [l^, revealing them as the outcome of nonlinear resonance (e.g. [sl, 
61] and references therein) among the various modes of the systems. Therefore, we argue 
that IFQPOs, which are about an order of magnitude smaller than HFQPOs/kHz QPOs 
(sometimes, however, IFQPOs are just about half the corresponding HFQPOs/kHz QPOs), 
also result from nonlinear resonance processes in accretion flows. 

If the accretion flow with radial and vertical epicyclic modes of oscillation is excited by 
an external mode, for the present model the disturbance mode due to the rotation/spin of 
ihe compact object, then the frequencies of original modes and disturbance mode are shown 
5|, |6| to have commensurable relations which may cause the modes to be strongly coupled 
and yield an internal resonance. This leads to the frequencies of newly formed modes having 
relations 

n m 

T^r + l^Vs = V, - — + —Us, 

when Ur + -v^ = Uu] ^ = (1) 

where and vi are the QPO frequencies, and Vz are respectively the radial and vertical 
epicyclic oscillation frequencies given by 



i/^ = ^^r2-4av^ + 3a2, (2) 
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where is the orbital frequency, A = — 2Mr + a^, r is the radial distance from the central 
compact object, Vg is the spin frequency of the compact object and m and n are integers. 
For a black hole, if the excitation and then corresponding coupling take place in such a way 
that m = n = 2, then from equation ([2]), Au = Vh — ^il^^s-, with Vz — Ur^Us- However, for a 
coupling with m = n = 1, Au = Uh — i^i ^ when z/^ — VrH^s/^. 



C. QPOs in GRS 1915+105 and its spin 

Figure [Iti shows that a pair of QPOs form for m = n = 2 due to the coupling among the 
modes yielding internal resonance. We know that GRS 1915+105 exhibits a pair of HFQPOs 
and its mass is expected to be in the range of 10 — ISMq. Hence, we choose M = IAMq in the 
above figure and the upper and lower HFQPOs are produced at = 168.5 Hz and ui = 113 
Hz respectively for a = 0.7, which are the observed HFQPOs for GRS 1915+105. Details 
are given in TABLE I. Note that each resonance produces two frequencies and hence IFQPO 
is absent in the above process. However, the source indeed exhibits an additional IFQPO 
with frequency z/j = 67 Hz. Hence, a successful model should reproduce this additional QPO 
frequency as well, for the same set of M and a which gives rise to the HFQPOs. We find that 
in fact the IFQPO with the observed z/j forms in our model for the coupling with m = n = 1 
of the fundamental modes, as shown in Fig. [T)d revealing internal resonance, along with 
z/;. However, if we choose M = ISMq, the extreme possible mass for GRS 1915+105, then 
z//i, ui, Pi are reproduced by our model for a = 0.77. See TABLE I for other details. 

Narayan and his collaborators, based on fitting the continuum spectra, argued jst that 
a > 0.98 for GRS 1915+105. Therefore, in our work, first we tried to confirm a large spin for 
this black hole. It, however, appears that a = 0.975 and 0.972, clearly < 0.98, to reproduce 
an IFQPO for M = 14Mq and 18Mq respectively for m = n = 1. Then for the pair of 
HFQPOs, our model predicts a = 0.95, again < 0.98, for the same m and at a choice of 



rather higher mass M = 18.4M( 



BE 



6|. Hence, not only a < 0.98, there is a clear mismatch 



between the sets of input parameters producing observed HFQPOs and IFQPO, questioning 



recently that the prediction of 



a high spin of GRS 1915+105. Indeed, it was shown 
spin based on diskoseismology and g-mode oscillation does not always tally with that from 
continuum fitting. Although for GRS 1915+105 the predictions from these two methods 
rather tally, for other sources there are problems. 
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TABLE I: M is given in units of Mq, Vh.i.i <^rc in Hz, locations of QPO {tqpo) and last stable 
circular orbit {rms) in accretion disk are in units of GM/(?. For GRS 1915+105, the observed 
Vh = 168, ui = 113, Vi = 67 and for XTE J1550-564, the observed Uh = 276, ui = 184, Ui = 92 
(controversial) . 

Sometimes, another QPO frequency ~ 41 Hz is reported for this source. If the set of 
frequencies 41 Hz and 67 Hz is considered to form a pair, similar to that having 113 Hz and 
168 Hz, then both the pairs should be reproduced by the same model. Assuming the first 
pair being the result of coupling with m = n = 1 and the second being with m = n = 2, 
they are reproduced for M — I6M0 and a — 0.62 in our model, which indicates the spin of 
GRS 1915+105 to be even slower. See TABLE I for other parameters. 
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D. QPOs in XTE J1550-564 and its spin 

The black hole source XTE J1550-564 also shows a pair of HFQPOs at 276 Hz and 
184 Hz along with an IFQPO at z/j = 92 Hz [16^. However, the IFQPO for this source is 
controversial. As the mass of the source (M = 9.1Mq) is also quite definite, we can test our 
model for this source in order to determine the spin of the black hole. In order to reproduce 
HFQPOs and IFQPO by the couplings with m = n = 2 and m = n = 1 respectively, as 
shown by Fig. [21 the spin of the source is found to be a = 0.715. Details are given in TABLE 
I. 

E. QPOs in fast spinning black holes 

Our model also predicts that a black hole with a larger spin is prone to exhibit single 
HFQPO based on a nonlinear resonance coupling with m = n = 1. Figure [3] shows a 
hypothetical case with a = 0.99 exhibiting, for m = n = 1, Uh ^ 330 Hz and a negative and 
hence unphysical ui. Similarly, single HFQPO forms for m = n = 2, except with a larger 
Ufi, which is very similar to that of a neutron star. Hence, a rapidly spinning black hole 
(and neutron star) seems to not favour the formation of a pair of HFQPOs. This is more 
unfavorable due to the fact that a fast spinning compact object presumably creates stronger 
disturbance in the surrounding disk, rendering only one QPO frequency dominant. 

F. QPOs in slow spinning black holes 

It is easy to check that for a low spinning black hole, the results with m = n = 1 and 
m = n = 2 (almost) coincide giving rise to low values of QPO frequency when Uh ~ i^i. 
Figure m shows that for a = 0.1, resonance occurs at a larger radius away from the black 
hole where merges to with — )■ 0. This is justified as, a smaller spin corresponds 
to a smaller strength of disturbance, which results in practically no nonlinear effects in the 
system, when the linear effect is also very weak. Hence, practically the coupling and hence 
the resonance takes place among the fundamental modes of the system — due to the radial 
and vertical epicyclic oscillations. This further implies that HFQPOs can not arise in low 
spinning black holes. 
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III. SPIN OF BLACK HOLES IN BLAZARS 

A. Classification of blazars 

Here we look at the issue of spin addressing the underlying luminosities. Blazars are 
the radio loud AGNs, which are divided into two classes: FSRQ and BL Lac. FSRQs are 
considered to be the subclass of FR II galaxies, whereas BL Lacs to be the subclass of FR I 



galaxies, based on the unification scenario 
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18|. Radio observations suggest that the jets 



of FR II galaxies are more collimated and more powerful than that of FR I galaxies. During 
its first year of observation, the Large Area Telescope (LAT) on board Fermi Gamma-ray 
Space Telescope (Fermi) detected around 600 blazars — almost an equal number of FSRQs 
(281) and BL Lacs (291) 
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20|. Based on their spectral energy distribution (SED), BL 



Lacs are further classified into three sub-categories: the low synchrotron peak (LSP), the 



intermediate synchrotron peak (ISP) and the high synchrotron peak (HSP) sources |21 |. 
With this significantly improved catalog, one can for the first time carry out a detailed 
comparison of 7-ray properties of FSRQs and BL Lacs and try to understand the underlying 
physics related to their intrinsic luminosities. Note that observationally FSRQs are on 
average about three orders of magnitude more luminous than BL Lacs. Also, BL Lacs 
exhibit a harder average photon spectrum than that of FSRQs and the luminosities of LSP 
sources are closer to the luminosities of FSRQs. From 7-ray spectral studies, it is observed 
that usually FSRQs show a break in the spectrum, while among BL Lacs, mostly LSPs 



show such a feature [2l|. Therefore, we consider FSRQs and LSPs as one group (FSRQ 
group) and ISPs and HSPs as another group (BL Lac group). In rest of the paper, we will 
use FSRQ group and FSRQs interchangeably and same for BL Lac group and BL Lacs. 
Figure [5] confirms that FSRQs exhibit on average three orders of magnitude higher observed 
luminosity and larger photon spectral index than those of BL Lacs. Moreover, FSRQs show 
strong emission lines, whereas BL Lacs do not. 

B. Properties of blazars 



Based on the isotropic luminosity and spectral index of these sources, earlier authors 22 1 
proposed that FSRQs accrete at or above a critical mass accretion rate (Mcric), whereas 
BL Lacs at a rate below Mcric, which results in their observed large luminosity differences. 
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However, as mentioned as well by others 23|, the actual cause of the FR I/FR II (broader 
classes of FSRQs/BL Lacs, as mentioned above) division is still not confirmed. There are 
two models. One argues that the interaction of jets with the ambient medium of different 
physical properties causes morphological differences between them (e.g. 2^). The other 
argues based on different intrinsic nuclear properties of accretion and jet formation processes 
(e.g., [25|) between them. 



Indeed, spectral energy distribution (SED) modeling suggests that jet emissions of FSRQs 
are dominated by external Compton (EC) processes. In this case, the relativistic particles in 
the jet upscatter external low-energy photons leading to the generation of 7-ray photons. The 
source of low-energy pho tons could be the underlying accretion disk, the broad-line region, 
the dusty torus etc. |26l-l35| . depending on how close the gamma-ray emitting region is to 
the black hole. On the other hand, SEDs of BL Lac sources are better fitted by synchrotron 
self-Compton (SSC) emissions, when the synchrotron photons of the jet are upscattered 
via inverse Compton processes by relativistic particles in the jet. Therefore, observed data 
suggest that the relativistic beaming effects in the two classes are not same and hence the 
intrinsic luminosity (which is the luminosity corrected for relativistic beaming), which is the 
property of the underlying accretion disk, should be treated as a more fundamental quantity 
than the observed luminosity in order to understand the origin of difference between the two 
source classes. However, if the source classes exhibit different accretion rates, intrinsic 
luminosities (i.e. that in the disk and hence at the base of the jet) also have to be different 
from each other. 

Mukhopadhyay and his collaborators j3] showed that the black hole spin plays a crucial 
role to control outflow and corresponding power. They considered a disk-outflow coupled 
region and by solving the hydrodynamic equations found that outflow is stronger for a faster 
rotating black hole than that for a slower one. Therefore, even with the same accretion rate, 
one can have a difference in luminosities due to the difference in the spin of central black 
hole. 



C. Difference between intrinsic and observed luminosities 

It is assumed that the blazar luminosity mostly comes from the jet which is believed to 
be made up of relativistic particles, either in a continuous flow or as discrete blobs. The 
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7-ray emission of jet arises due to either SSC and/or EC processes, as mentioned above. 
The observed beamed luminosity {Lots) and the jet frame luminosity {Lunbeamed) are related 

by laei 

Lobs Lunbeamed ^ ; (3) 

where m = 2 for a continuous jet and m = 3 for a discrete jet, n = for SSC process 
and n = 2a^ + 1 for EC process, when is the energy spectral index. One can obtain 
the intrinsic luminosity Lunbeamed of the source from equation ([3]) by knowing Lohs and 5. 
We use an average value of 5 for blazars and individual spectral indices while calculating 
unbeamed/intrinsic luminosities of FSRQs and BL Lacs. In order to obtain Lunbeamed, we 
consider a continuous jet for the BL Lac group having SSC emission. However, for the 
FSRQ group, we consider a combination of SSC and EC emission models in the framework 
of continuous jet. Since the fraction of SSC and EC combination is not well constrained, we 
consider different scenarios: (a) 50% SSC and 50% EC contributions, (b) 25% SSC and 75% 
EC contributions, (c) 10% SSC and 90% EC contributions, and (d) 100% EC contribution. 



D. Fixing 6 in order to determine intrinsic luminosities 

A very important factor to determine the intrinsic luminosity of a blazar is 6. There 
are mainly two approaches to determine 6, actually Lorentz factor F which is related to 
6: (i) SED modeling, (ii) radio Very Large Baseline Interferometry (VLBI) measurement. 
Considering the opacity arguments, it is generally believed that the 7-ray emitting region is 
located closer to the central black hole. Hence, SED modeling beyond radio may be more 
appropriate, since it more accurately represents the region of 7-ray emission. Unfortunately, 
in most cases SED data gathered across wavelengths are not simultaneous. Further, SED 
modeling involves a large number of model parameters and, hence, typically results in large 
uncertainties. Another important parameter, jet to line-of-sight angle {9 jet)-, is not measured, 
rather a nominal value of 9 jet is typically considered for this modeling. Moreover, for most 
of the cases, the goodness of fit is not examined and, hence, the parameter values are not 
optimized. Based on SED modeling of 85 blazars from the first 3 months' Fermi observation, 

n 

it was found [37[ that the values of F for blazars are within the range of 10 — 15. 

The values of F and 5 can also be determined from VLBI observations based on direct 
observation of the brightness temperature of the source {Tb^obs)- Comparing T^^obs with the 
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intrinsic brightness temperature of the source (Tb j„t), one can determine S 38|, l39|. From 



the observed apparent velocity (vap) and 6, Ojet and V can be derived. However, in this 
technique Tb^jji* is o 



earher works 
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"ten assumed to be the equipartition temperature. It was considered in 
39| that Tij^int is same for all sources. This assumption makes the estimated 
values of F and 5 less accurate. 

The values of F and 6 could also be determined \^\ from VLBI observations adopting yet 
another technique. In this technique, 5 is derived by comparing the timescale of decline in 
flux density with the light travel time across the emitting region. From the knowledge of 5 
and Vapi the previous authors j^oj derived 9 jet and F. This approach yields more appropriate 
values in comparison to other alternate methods. In order to determine the average 5 
and F from this approach |4o|, we assume that the intrinsic 5 (also F) distribution and 
the corresponding errors are Gaussian (following the methodology used in |41j] and 42|] to 
determine the average spectral index of blazars). The average 5 for FSRQs is then obtained 
to be 23.1 ± 8.9 and that for BL Lacs to be 15.3 ± 5.5 and the respective values of F to 
be 17.2 ± 5.3 and 12.5 ± 3.5. For the present work, we choose 5 and F obtained from this 
approach. Considering the large errors, we adopt the average values oi 5 = 20.6 ± 8.4 and 
F = 15.1 ± 4.6 for both the source classes. 

There is, however, a concern regarding the difference in location of 7-ray emission versus 
the location of radio emission. However, the independent estimates of F values derived 
from the SED modeling js^] lie in the range 10 — 15, which is consistent with the mean 
value adopted here. This probably indicates that the parameters 5 and F are not varying 
significantly between the two regions. We use the average value of 5 = 20.6 to estimate the 
unbeamed luminosities. 



E. Evaluating intrinsic luminosity and its dependence on the spin of black hole 

With the above 5, using equation (|3]), interestingly the unbeamed/intrinsic luminosity 
difference between FSRQs and BL Lacs is shown in Fig. [6] to decrease significantly in com- 
parison to the observed values. There is, however, a small difference in intrinsic luminosities 
between BL Lacs and FSRQs (with FSRQs on average having higher luminosity). Neverthe- 
less, the small difference in intrinsic luminosities suggests that it is difficult to accommodate 
models involving vastly different accretion rates in these two systems. Any change in ac- 
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cretion rate would affect the properties of the flows in the accretion disk itself, before it is 
launched and beamed. We propose that the intrinsic luminosity mismatch arises from the 
difference in the spin of the central black holes. Indeed one expects some connection between 

nn Q 

the jet and the central black hole [43|, |4^. In an earlier paper [7[, based on a disk-outflow 



coupled model, the fundamental properties of central black hole were shown to be linked 
with accretion and outflow, which in turn could be related to the properties of FSRQs/BL 
Lacs. In that work, the total mechanical outflow luminosity was found to increase with the 
increase of spin of the central black hole. This implies, for a faster spinning black hole, 
matter will flow out faster, making the outflow stron ger with higher flow density. 

For a fixed accretion rate, it is known (see e.g. for the latest solutions) that 

with the increasing spin of black hole, the optically thick part of the accretion disk advances 
towards the black hole. Therefore, the supply of soft photons to the jet emitting region 
increases with the increasing spin of black hole. On the other hand, as shown earlier 
with the increase of spin, density of the disk-outflow coupled region at a particular accretion 
rate increases. Therefore, with the increase of spin, denser matter having larger supply of 
soft photons will interact with radiation more rapidly, which results in losing more energy. 
Hence matter in the jet and then corresponding radiation is expected to exhibit a steeper 
spectrum. This agrees with the observed steeper 7-ray spectrum of FSRQs than that of BL 
Lac objects. It was found earlier ^|, based on the host galaxy optical luminosity — radio 
luminosity plane dividing FR I and FR II radio galaxies, that radio power is proportional 
to the optical luminosity of the host galaxy. However, as discussed in §111. D that basic 
properties in radio and 7-ray emitting regions do not vary significantly. Hence, all the 
above properties also explain the existence of strong ionizing continuum in FSRQs, which 
are expected to harbor faster spinning black holes (as explained further in §111. F below) and 
therefore exhibit denser accretion flow, revealing strong emission lines in them. 



F. Spin difference between black holes in FSRQs and BL Lacs 

It is also very important to note that although the unbeamed luminosities of FSRQs 
and BL Lacs are much smaller than their observed values (and the difference of unbeamed 
luminosities between them decreases significantly, as described above), FSRQs still exhibit 
higher luminosities than BL Lacs. This implies that FSRQs harbor fast spinning black holes 
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than BL Lacs [Tj]. Note that it is unreasonable to beheve that there is a 100% EC contribution 
in observed luminosities of FSRQ group, which will make their intrinsic luminosities and 
hence spin to be smaller than that of BL Lacs — that goes against their spectral behavior 
as mentioned above. Therefore, we consider the bound on the contribution from EC process 
to be <95%. 



If one considers that the magnetic field plays a major role in collimating the jets 44j. |48| . 



it can result in more twisting of magnetic field lines and, hence, a more coUimated jet for 
a faster rotating black hole, lending further support for enhanced spin in FSRQ systems. 
By numerical simulations, it has been already investigated [S] that whether the spin of 
central black hole can play a major role to explain radio loud/quiet dichotomy or not. For 
a geometrically thick disk, as is also considered in p], simulations show that the choice of 
a spin value a~0.15 of a black hole for radio quiet AGNs can explain their difference in 
luminosities from radio loud AGNs having a maximally spinning black hole. 



G. Predicting spin for black holes in FSRQs and BL Lacs 

Now recalling the relationship connecting the total mechanical power of outflow (Pj) with 
the spin of black hole given by equation (20) and shown in Fig. 7(a) of [7|, we derive a simple 
relationship between Pj and a, as shown by Figure [Tl^a), based on fitting the power given 
by the above mentioned equation in jr^ by a third order polynomial of a as 

where A = 2.87±0.26, B = -4.08±0.40, C = 2.88±0.17 and D = 41.53±0.02. Since the 
earlier work [7] did not consider any beaming effect (F = 1) while calculating Pj, which is 
calculated in the frame of disk-outflow coupled system, we assume in the first approximation 
that Pj is proportional to the intrinsic luminosity such that 



Pj,FSRQ _ Lunbeamed,FSRQ 



(5) 



p r 

j,BLLac J-'unbeamed,BLLac 

but not the observed luminosity, due to the inadequate knowledge of the underlying radiation 
hydrodynamics therein. 

Now, for a particular spin of the black hole in BL Lac, we calculate Pj^bllgc using equa- 
tion (jl]). Since the ratio of Pj of the two systems is equal to the ratio of their intrinsic 
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luminosities, the corresponding power of FSRQ (Pj^psBQ) can be estimated by knowing the 
average intrinsic luminosities of the two classes. Therefore, the corresponding average spin 
of the black hole in FSRQ is derived from equation (jlj). We compute the average spin of 
FSRQs' black hole for a range of the average spin of BL Lacs' black hole, as shown in Fig- 
ure [71(b). FSRQs and BL Lacs, being radio loud AGNs, are expected to be harboring fast 
rotating black holes than the radio quiet AGNs 8| . Hence, assuming that the spin of black 
holes in radio loud AGNs (and hence of BL Lacs, which are argued above to be harboring 
slower black holes among blazars) should not be less than 0.5, the minimum possible spin for 
the black hole in FSRQs is estimated to be ~0.9, for a case where the SEDs of FSRQs are 
assumed to have 75% EC contribution (and ~0.68 assuming 90% EC contribution). Con- 
versely, assuming FSRQs to be harboring maximally spinning black holes, the maximum 
possible spin for the black hole in BL Lacs is ~0.75, considering 75% EC contribution (and 
~0.93 for 90% EC contribution). 



IV. DISCUSSION AND CONCLUSION 



By simultaneously modeling three classes of QPOs of the black hole in XRBs — a pair of 
HFQPOs and an IFQPO, we predict the spin of black holes. While the model can predict 
spin by fitting either a pair of HFQPOs or an IFQPO (for black holes exhibiting only one 
of the classes), the result will be robust for the black holes showing both classes of QPO. 
The sources GRS 1915+105 and XTE J1550-564 are the ones exhibing both HFQPOs and 
IFQPO and we predict their spins. Our model predicts that GRS 1915+105 could be fast 
spinning (a = 0.95 — 0.972), only if its mass M>18Mq. But this set of mass and spin could 
not predict the existence of HFQPOs along with IFQPO. For the existence of both HFQPOs 
and IFQPO, our model argues for a = 0.77, for the choice of M = ISM©. For M = 14M0, 
a is obtained to be 0.7, which is a more robust result. 

The present model also argues that a fast spinning black hole cannot exhibit a pair of 
HFQPOs; it reveals only one HFQPO. Moreover, a slow spinning black hole cannot exhibit 
any HFQPO (and presumably IFQPO). Hence, the pair of HFQPOs in a 2 : 3 ratio seems 
to be the property of a black hole spinning neither very fast nor slow. Same could be true 
for neutron stars exhibiting kHz QPOs. 

In principle, above ideas could be applied for SuBHs in oder to determine their spin. 
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Unfortunately, due to their very low expected frequencies, it is very difficult to probe their 
QPOs. Hence, in order to predict the spin of SuBHs, we model the sources, for the present 
purpose blazars, differently. 

We propose that there is a difference in the spin of the central black holes in the two 
classes of radio loud AGNs, namely, FSRQs and BL Lacs. The spin difference between these 
two classes of radio loud AGNs is predicted based on their difference in the mechanical 
outflow powers (and hence, the unbeamed/intrinsic luminosities). We predict that FSRQs 
are expected to harbor a faster spinning black hole than BL Lacs. By the same argument, 
the black hole spin in LSPs is expected to be higher than that in ISPs, while that in HSPs 
to be the slowest. However, both the source classes must be harboring spinning black holes. 

Unfortunately, due to lack of proper data required for this modeling, we could not predict 
the spin of individual blazars. We rather attempt to predict the average spin of FSRQs and 
BL Lacs from the large sets of respective data and show that general relativistic effects of 
rotating black hole play important role in these systems. 

Note that recently a new type of jets was observed in the sources of SuBHs due to tidal 
disruption of stars which are supported ouly .,v the existence of rap.dly spiuu.ug black 
holes at the center [50 1. 

More detailed observations of these systems in future will provide a potential laboratory 
to investigate the underlying general relativistic processes of these sources. 

Therefore, by modeling observed data from XRBs and blazars based on the principle 
of general relativity, we predict the spin of the black holes in the respective sources to be 
nonzero. This provides a natural proof for the existence of Kerr metric. 

Acknowledgments: This work was partly supported by an ISRO grant 
ISRO/RES/2/367/10-11. Thanks are due to Upasana Das for carefully reading the 
manuscript. 
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FIG. 1: Representative case of GRS 1915+105: (a) Variation of frequencies of the higher (dotted 
hne) and lower (dashed hne) HFQPOs, their difference (long-dashed line) and difference between 
vertical and radial epicyclic frequencies (dot-dashed line) in Hz as functions of radial coordinate 
in units of GM/c^ of the accretion disk for m = n = 2. (b) Variation of frequencies of the 
lower (dotted line) of HFQPOs and IFQPO (dashed line), their difference (long-dashed line) and 
difference between vertical and radial epicyclic frequencies (dot-dashed line) in Hz as functions of 
radial coordinate in units of GM/(? of the accretion disk for m = n = 1. The vertical solid line 
pinpoints the location of resonance. Other parameters are M = 14Mq, a = 0.7; see TABLE I for 
other details. 
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FIG. 2: Representative case of XTE J1550-564: (a) Variation of frequencies of the higher (dotted 
Une) and lower (dashed hne) HFQPOs, their difference (long-dashed line) and difference between 
vertical and radial epicyclic frequencies (dot-dashed line) in Hz as functions of radial coordinate 
in units of GMjc? of the accretion disk for m = n = 2. (b) Variation of frequencies of the 
lower (dotted line) of HFQPOs and IFQPO (dashed line), their difference (long-dashed line) and 
difference between vertical and radial epicyclic frequencies (dot-dashed line) in Hz as functions of 
radial coordinate in units of GM/(? of the accretion disk for m = n = 1. The vertical solid line 
pinpoints the location of resonance. Other parameters are M = %.1Mq, a = 0.715; see TABLE 1 
for other details. 
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FIG. 3: Representative case of a high spinning black hole with a = 0.99: (a) Variation of frequen- 
cies of the higher (dotted line) of HFQPOs, the difference between higher and lower (unphysical) 
frequencies of HFQPOs (long-dashed line) and difference between vertical and radial epicyclic fre- 
quencies (dot-dashed line) in Hz as functions of radial coordinate in units of GM/(? of the accretion 
disk for m = n = 2. (b) Same as in (a), but for m = n = 1. The vertical solid line pinpoints the 
location of resonance. M = IbM^; see TABLE I for other details. 
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FIG. 4: Representative case of a low spinning black hole with a = 0.1: (a) Variation of frequencies 
of the higher (dotted line) and lower (dashed line) HFQPOs, their difference (long-dashed line) and 
difference between vertical and radial epicyclic frequencies (dot-dashed line) in Hz as functions of 
radial coordinate in units of GM/(? of the accretion disk for m = n = 2. (b) Same as in (a), but 
for m = n = 1. The vertical solid line pinpoints the location of resonance. M = ^Mq; see TABLE 
I for other details. 
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FIG. 5: Distribution of observed luminosity and spectral index for FSRQs and BL Lacs. Red 
points represent FSRQs, green, blue and violet points respectively represent LSPs, ISPs and HSPs. 
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FIG. 6: Same as in Fig. [5l except intrinsic/unbeamed luminosities are plotted in the horizontal 
axis. Contributions from 50% SSC and 50% EC are chosen for FSRQs. 
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FIG. 7: (a) Total mechanical power of the outflow Pj (points) as a function of spin (from 
when solid line represents the fit with the function given by equation ([4]). (b) The values of FSRQ 
spin for a range of BL Lac spin. We take out the beaming effect considering SSC emission in BL 
Lacs and a combination of SSC and EC emissions in FSRQs. 
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